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Abstract: Motivated by the observed cosmic matter distribution, we present the following 

"-^ ' conjecture: due to the formation of voids and opaque structures, the average matter density 

I ■ on the path of the light from the well-observed objects changes from ^1^/ ~ 1 in the 

;_( . homogeneous early universe to ^m — in the clumpy late universe, so that the average 

c/3 , expansion rate increases along our line of sight from EdS expansion Ht ~ 2/3 at high 

C^ . 

redshifts to free expansion Ht ~ 1 at low redshifts. To calculate the modified observable 



distance-redshift relations, we introduce a generalized Dyer-Roeder method that allows for 



K*" ' two crucial physical properties of the universe: inhomogeneities in the expansion rate and 

i^ ' the growth of the nonlinear structures. By treating the transition redshift to the void- 

CN ' dominated era as a free parameter, we find a phenomenological fit to the observations from 

the CMB anisotropy, the position of the baryon oscillation peak, the magnitude-redshift 

relations of type la supernovae, the local Hubble fiow and the nucleosynthesis, resulting in 

C — ! a concordant model of the universe with 90% dark matter, 10% baryons, no dark energy, 

^^ ' 15 Gyr as the age of the universe and a natural value for the transition redshift zq = 0.35. 

Unlike a large local void, the model respects the cosmological principle, further offering 

^ ' an explanation for the late onset of the perceived acceleration as a consequence of the 

C^ ■ forming nonlinear structures. Additional tests, such as quantitative predictions for angular 

deviations due to an anisotropic void distribution and a theoretical derivation of the model, 

can vindicate or falsify the interpretation that light propagation in voids is responsible for 

the perceived acceleration. 
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1. Introduction 

The widely adopted framework of modern cosmology is the spatially homogeneous and 
isotropic Friedman-Robertson- Walker spacetime, with the growth of structure described as 
linear perturbations evolving on the smooth FRW background. Accordingly, the standard 
model of cosmology is built on the assumption that the effects of the evident nonlinear 
inhomogeneities on the detectable light average out^ over cosmological distances. This as- 
sumption - although perhaps in concordance with Newtonian intuition ~ lacks a convincing 
demonstration within general relativity, not to mention a mathematically rigorous proof. 



^Apart from gravitational lenses that are occasionally taken into account as additional corrections. 



In fact, the assumption was questioned already in the 60's by Zel'dovich [1], Feynman^, 
Bertotti [2], Gunn [3], Kantowski [4] and in the 70's by Dyer & Roeder [5,6], who sug- 
gested that the dumpiness of matter in the real universe has consequences on the observable 
distance-redshift relations that the simplified FRW description fails to capture. 

At the time of this pioneering research, the cosmological observations were too inaccu- 
rate to distinguish the predictions of the proposed inhomogeneous models from the simple 
and elegant FRW solutions. Consequently, the efforts towards a more thorough description 
of the effects of the nonlinear structures on the distance-redshift relations never became 
part of the mainstream research. Although the structure formation itself has since then 
developed into a major part of modern cosmology, the emphasis is usually given either to 
the general relativistic but perturbative description of structures on large scales, or to the 
Newtonian description of the small scale structure e.g. in the context of galaxy formation. 

At present, the unparalleled precision of the recent cosmological observations has raised 
the debate about the effect of the nonlinear structures on the distance-redshift relations 
once again topical. Indeed, the observations indicate the existence of inhomogeneities that 
seem to be on large enough scales to have cosmological significance but too strong to be 
described within the linear theory — in particular the observed voids and the clustering of 
matter as filaments of opaque galaxies between the voids [7-12]. The increase in the preci- 
sion of the data naturally demands a corresponding improvement in the theoretical model; 
otherwise there is an increased risk of drawing incorrect conclusions from the observations. 

In fact, the discovery that the standard FRW cosmology needs a fine-tuned dark energy 
component [13-15] in order to account for observations [16-19], suggests that the precision 
of the data may already have surpassed the precision of the employed theoretical model. 
This is also backed up by a mysterious coincidence, pointed out by Schwarz [20] and 
Wetterich [21] and elaborated by Rasanen [22-25]: the perceived increase in the cosmic 
expansion, ascribed to dark energy in the standard cosmology, happens to take place during 
the same era when nonlinear structures start to form at cosmologically significant scales. 

Whereas the fine-tuning of dark energy only indicates that something fundamental 
may be lacking from the standard picture, the coincidence problem suggests a causal con- 
nection between cosmic acceleration and structure formation. Perhaps the most natural 
explanation would be, that the standard FRW description breaks down when the nonlin- 
earities become dominant in the late universe and the need for a tiny cosmological constant 
is a manifestation of this breakdown, not evidence for dark energy [22]. In recent years, 
more and more authors have addressed this issue and substantially increased the general 
understanding of the subject, but crucial open questions still remain [26]. In this work, 
we address the issue by considering a phenomenological model that aims to: 1) clarify the 
physical interpretation for how the structure formation can actually mimic dark energy, 
and 2) provide a quantitative fit for the observations. 

The paper is organized as follows. In Sects. 2.1-2^, we review three physically inde- 



pendent mechanisms for how the cosmic inhomogeneities can induce accelerated expansion 
along our line of sight: global acceleration due to spatial variations in the expansion rate, 

^Apparently, R. P. Feynman was also one of the pioneers of this subject, as [3,4] cite a colloquium by 
Feynman in 1964. 



faster local expansion rate due to a large local void, and biased light propagation through 
voids that expand faster than the average; possibilities to describe the total effect of the in- 
homogeneities have been discussed in Sect. ^^. After reviewing the original Dyer-Roeder 



method in Sect. 2.5, and motivating the need to modify it, we introduce a generalized 



Dyer-Roeder method in Sect. 3.2 to allow for the total effect of the inhomogeneities on the 
observable distance-redshift relations. The distance-redshift relations are then calculated 
for an inhomogeneous model, based on the conjecture that the average matter density on 
the path of the light from the observed objects changes from Qm — 1 in the homogeneous 
early universe to fijv/ — in the clumpy late universe. The physical foundations of the 
conjecture are summarized in Sect. |3.1| . In Sect. ^^, we construct a simple description 
of the structure formation in terms of a function a{z) that measures the deviation of the 
matter density along our line of sight from the FRW matter density. A comparison of the 
model with cosmological observations is performed in Sect. |3.4|, and a discussion of the 



results is given in Sect. 3.5. Finally, Sect. H contains our conclusions. 



2. Inhomogeneous generalizations of the FRW cosmology 

Formally, the cosmic inhomogeneities pose the following problem: we have a physical sys- 
tem, the universe, with > 10®'^ degrees of freedom but we want a model with < 10 degrees 
of freedom, represented by the cosmological parameters. For a mathematical description 
of the universe, we thus need a coarse graining map: 10®'' dof ^^ 10 dof , that preserves the 
relevant physical structure of the spacetime but removes its mathematical complexity. As 
our cosmological knowledge is firmly based on observing light, the relevant structure is the 
correct optical properties of the universe. That is, we want the observable distance-redshift 
relations to be unaltered by the coarse graining. 

A common critique against acceleration from inhomogeneities is that linearly perturbed 
FRW solutions seem to describe well even highly nonlinear density contrasts, Sp/p ~ 10'^'', 
while still having insignificant backreaction [27]. This is because the relevant perturbed 
gravitational potential can satisfy <1> ^ 1 even for the high density contrasts dp/p ^ 1. 
However, this argumentation overlooks the crucial part of realistic nonlinear structures, 
namely inhomogeneities in the expansion rate [23]. Even a relatively small spatial varia- 
tion in the expansion rate, 69/6 ~ 0(0.1), can invalidate the linear perturbation theory by 
rendering the uniformly expanding FRW background inapplicable. Sometimes the assump- 
tion 69/6 ^ 1 is hidden in the requirement of small peculiar velocities [28], which appears 
to restrict the spatial variations in the expansion rate to be small as well [29]. As we 
next discuss, O(O.l) spatial variations in the expansion rate can indeed induce accelerated 
expansion via various mechanisms. 

2.1 Backreaction due to inhomogeneous expansion 

A central question in constructing the coarse grained description of the universe is how 
to average relativistic gravitational systems within Einstein's theory. As pointed out by 
Shirokov and Fisher [30] and later made more popular in the observational cosmology 
programme by Ellis and collaborators [31-34], it seems physically more correct to first 



calculate the Einstein field G{g) for the exact metric g and only then average {G{g)), than 
to calculate the Einstein field for an averaged metric G{{g)) as in the standard FRW cos- 
mology. The reason is that the Einstein field is more closely related to physical quantities 
whereas the metric corresponds to a gravitational tensor potential, whose derivatives de- 
termine the physics. In general relativity, the field G depends nonlinearly on the metric 
g, so its evaluation does not commute with averaging and the non-vanishing difference, 
{G{g)) — G{{g)) ^ 0, gives rise to what is known as the nonlinear backreaction; see [22-24]. 
Hence the issue is not only a conceptual one, but in general, the two approaches yield 
identical predictions only in the absence of nonlinear inhomogeneities. 

Even if performed in the above explained order, the averaging is not free of problems. 
Firstly, the only practicable way to average tensors seems to be the manifestly coordinate- 
dependent averaging of components: (G) = {G^y). As we would like the coordinate- 
invariance of general relativity to hold also for the coarse-grained or averaged quantities, 
only the inherently invariant rank tensors or scalars appear to have well-defined averages. 
Secondly, the desire in cosmology to divide the universe into temporal and spatial parts 
leads to the issue of foliation dependence: in order to take spatial averages, we must 
artificially break the symmetry between time and space inherent in general relativity. 

In 1999, Buchert proposed a way to circumvent, or at least alleviate, these problems 
[35]. By using the fact that a part of the Einstein tensor equation can be decomposed 
into scalar equations, he derived a set of spatially averaged equations, now known as 
the Buchert equations, describing the averaged dynamics of a general irrotational dust 
universe. Moreover, he foliated the spacetime by using the proper time of the dust as a 
global time coordinate; a physically well justified choice in a dust universe, also supported 
by a comparison with exact observables performed in Ref. [36]. 

To see how the backreaction can affect the observations, consider the Buchert acceler- 
ation equation [35]: 

3- = -47rG(p) + Q , (2.1) 

a 

where a{t) is a generalized scale factor, p is the matter density and the difference between 
Eq. (|2.1| ) and its homogeneous FRW counterpart is known as the dynamical backreaction 

Q^lm-{e)')-2{a'), (2.2) 

where o"^ is the shear scalar and 9 is the expansion scalar. The backreaction ( |2.2| ) in the 
Buchert equations explicitly demonstrate that averages of inhomogeneous quantities do not 
evolve in time like the corresponding homogeneous quantities. 

From Eq. ( |2.1j ) it is obvious that, regardless of decelerating locally everywhere, with 
large enough backreaction Q > 4ttG{p), the average expansion can accelerate without an 
exotic fluid with negative pressure or a modification of gravity. Physically, the global ac- 
celeration is possible because the volume can become dominated by fast-expanding regions 
[24]. This is realized with large enough variance of the expansion rate, if the counter- 
balancing average shear is not too large. The variance becomes large when contracting 
{9 < 0) and expanding (9 > 0) regions coexist. This is exactly what one would expect in 



the late universe with structures forming via gravitational collapse [37,38], so it has been 
conjectured that the average acceleration could explain the cosmological observations [23]. 

To clarify the essential point, consider two disconnected decelerating regions U and 
O that have initially the same volumes Vi. Let the region U represent an underdensity 
that initially expands at the rate 9i > and let O represent an overdensity that initially 
contracts at the rate —6i < 0. Clearly, the volume-averaged expansion rate is initially zero: 
{d)i = [ViOi + Vi{—9i))/2Vi = 0. Later, the region U has expanded to take up a volume 
Vu > Vi which is still expanding, 9u > 0, though slower than initially, 6u < 9i, whereas 
the region O has shrunken to a small virialized structure with volume Vq <^ Vu which is 
essentially static, Oq — 0. In the process, the volume-averaged expansion rate has become 
positive: {9) = {Vu9u + Vo9o)/(yu + Vq) ^ 9u > 0, confirming that the global expansion 
has accelerated regardless of decelerating locally at each point. By taking the junction 
conditions between the different regions into account, the emergence of the acceleration 
has been verified within exact solutions of general relativity [38, 39]. 

The current studies suggest that while the dynamical backreaction (^]^) can be sig- 
nificant, perhaps increasing the average expansion by a factor ~ 1.3, it may not be large 
enough to account for the required total increase factor ~ 1.5 alone [40-45]; however, it has 
been argued that the approximations in the current calculations may underestimate the 
backreaction [24]. In Sect. ^^ we propose that, due to the formation of opaque structures, 
the volume-averages do not necessarily capture the total effect of the inhomogeneities but 
the cosmic expansion may increase more rapidly along our line of sight to the observed 
objects than the global volume- average. 

2.2 A large local void 

When looking at the large-scale structure of the universe, perhaps the most eye-catching 
feature is that the major part of the universe appears to be taken up by voids of size 
0(10).. .0(100) Mpc [7-11], nearly empty regions expanding faster than the whole universe 
on average. Even larger voids could exist: it has recently been suggested that the correlation 
between the cold spot in the CMB and the deficit of radiogalaxies in the same direction 
might be evidence for a void of size ~ 300 Mpc [12]; though see [46]. It could also be 
possible that we would happen to live inside such a void [47-50], also known as Hubble 
Bubble in this context. 

Clearly, if we were sitting in a special spatial location, such as inside a large local void, 
the global spatial averages may fail to describe our cosmological observations. In that 
case, an option^ is to give up the notion of a global scale factor and study an inhomoge- 
neous metric; a particularly useful metric for describing voids is the spherically symmetric 
Lemaitre-Tolman-Bondi metric, an exact solution of the Einstein equations that Lemaitre 
discovered in 1933 [51]; for a review, see [52]. During the recent years, the LTB solution 
has been actively exploited in studying voids and the effects of inhomogeneities on the 
cosmological observations [47-49,53-63]. 

A key notion in understanding the physical basis of the inhomogeneity induced per- 
ception of acceleration is that the cosmological observations are made along our past light 
^Another option could be to consider a different averaging scale for each redshift [36]. 



cone. A way to illustrate this is to consider the directional derivative along the light cone 

dz dzdx^ dzdt^ ° V^^ dt) ' ^'^ 

where the approximation in the last step is more accurate for small distances, r <C Hq , 
but the sign is correct even for larger r. Eq. ( |2.3| ) tells us that negative radial variation 
roughly corresponds to positive time variation and that their relative contributions to the 
total variation w.r.t. the observable redshift cannot be distinguished. This is natural since 
by looking at a source, we simultaneously look into the past and spatially further. 

Within the LTB model, the argument can be made more precise by comparing the 
expression of the luminosity distance in a matter dominated local void [36] 



di™(z) = i7o-i(0) 






where ' = d/dr, with its expression in the homogeneous and flat ACDM model: 



d^--(z) = H, 



^+ [l + ^-nA z^ + [A-n^ + ^nl] z^ + oiz' 



(2.4) 



(2.5) 



The comparison between Eqs. (2^) and ( ^.51 ) confirms that negative radial variation of the 
expansion rate, Hq{0) < 0, is needed to mimic positive cosmological constant, f^A > 0, or 
to induce accelerated expansion along our line of sight. In other words, for the expansion 
rate to increase towards us along the past light cone, the LTB expansion function Ho^r) 
must decrease as r grows. This is exactly what an observer inside a void would see, and 
has been shown to account for various cosmological observations without local acceleration 
of the expansion [55, 56, 58, 60-62]. 

However, the explanation of the perceived cosmic acceleration as an effect induced 
by a local void does not come for free. Firstly, to respect the cosmological principle, we 
should not sit in the center of a spherically symmetric void, like fitting a local void for 
the observations seem to require; e.g. for the LTB model to remain consistent with the 
observed CMB dipole without a fine-tuned peculiar velocity, we should reside within a few 
percent of the size of the void from its center [64]. In addition, even if the void were not 
perfectly spherical and even if we were not near the center of the void, the cosmological 
principle would still be, at least weakly, violated: there are less galaxies in a void and if 
born in a random galaxy, it is much more probably located in an environment with lots 
of galaxies than in a relatively empty void. Another issue is that the required size of the 
local void seems to be a few times bigger than the typical size of the observed large voids. 
Moreover, the existence of such a void has not been observationally confirmed, although, 
due to most of the matter in the universe being dark, it may not be easy to disprove it 
either. In any case, without a natural explanation for these issues, the scenario with a large 
local void - if taken literally - does not appear likely. Instead, if considered to describe 
observations averaged over the celestial sphere, the LTB solution may still have relevance 
as an effective description for light propagation through several smaller voids. 



2.3 Selection effects due to opaque structures 

Even if we were not located inside a large local void, the observed cosmic network of 
smaller voids can still significantly affect the observable distance-redshift relations. Firstly, 



as discussed in Sect. 2.1, the spatial variations in the expansion rate can increase the 
volume-averaged global expansion rate. In addition, a complementary aspect we want to 
point out in this work is that the opacity of the cosmic structures can bias the observations, 
such that the light we see has traveled mostly in freely expanding voids at late times, and 
thereby affect the observed distance-redshift relations, independently from the dynamical 
backreaction. Although it may be difficult to distinguish the observational consequences 
of the selection effects and the backreaction, they are in principle distinct mechanisms: 
one could have large opaque structures that induce selection effects but give negligible 
backreaction or, on the other hand, perfectly transparent but highly nonlinear structures 
that produce backreaction but no selection effects. Naturally, in the real universe, both 
mechanisms may have significant effects on the observations. Mathematically, the difference 
between the purely dynamical backreaction and the biased light propagation in voids can 
be seen as the difference between the volume- averaged expansion rate and the average 
expansion rate along our line of sight. 

The selection effects in the observations of cosmic light can be justified by the fact that 
the emptier a region in space, the easier it is to see an object through it. The bias can be 
further categorized roughly as two different cases: 1) an opaque structure on the foreground 
simply screens the further object behind it, or 2) the further object is seen but, due to 
foreground contamination, not clearly enough that the object would contribute to the 
high quality data. The first case is perhaps relatively rare, since only the nearest galaxies, 
including the Milky Way, cover a notable fraction of the celestial sphere. Instead, the second 
case may be more relevant as, due to their long distance, the cosmological objects tend to be 
dim, which means that they can be harder to detect if there are bright foreground objects 
nearby on the sky. In addition to the brightness of the foreground, the contamination can 
arise from gravitational lensing. An important implication of the contamination is that 
the effective area screened by the foreground objects may be significantly larger than their 
physical cross-section. As an example, large type la supernova surveys typically cover only 
~ 10~^ part of the sky [65]. It thus appears conceivable that the well-observed supernovae 
may consist more of objects with less than average intervening material on their foreground. 

In addition to the kind of bias discussed above, there may be other kind of selection 
effects, arising from the fact that the dimmer the object the harder it is to detect. It appears 
non-trivial how this would affect the observed distance-redshift relations, since on one hand 
it could favor objects with higher luminosity but on the other hand disfavor objects with 
high redshift. In any case, for the statistically most relevant brighter supernovae at lower 
redshifts, we would expect this effect to be less significant. Furthermore, if the type la 
supernovae would be perfect standard candles, the bias towards higher intrinsic luminosity 
should not exist at all. In addition to selection effects, there could also be other systematic 
effects related to observational techniques [66]. Here we simply assume that such effects 
do not play a significant role. 



Unfortunately, due to the complexity of the cosmological structures, it appears ex- 
tremely difficult to derive quantitative estimates for the significance of the selection effects. 
A possible way to obtain such estimates could be a computer simulation where one studies 
the difference of the average expansion rate along our line of sight relative to the volume- 
averaged global expansion rate. Altogether, we can currently neither vindicate nor rule 
out the importance of cosmic selection effects from opaque structures. 

2.4 Total effect of the inhomogeneities 



As we have discussed in Sects. |2.lH 2.3|, it appears very difficult to derive from theory a 
realistic inhomogeneous model that would take into account all the relevant effects of the 
inhomogeneities. Instead of dealing with these difficulties, we employ a phenomenolog- 
ical approach based on physically justified assumptions in Sect. ^ The model contains 
a transition redshift to the void-dominated era as a new free parameter and the applied 
mathematical framework is provided by a generalized Dyer-Roeder method, introduced in 



Sect. 3.2. In this approach, the observable or optical distance-redshift relations are deter- 
mined by the average properties of the space along our line of sights to the objects, instead 
of global volume- averages. The aim of the model is to provide only a simple description of 
what the overall effect of the cosmic inhomogeneities on the observed distance-redshift re- 
lations may be. Although based on assumptions, in the end, one can use the naturalness of 
the required value for the new parameter to judge the viability of the model, e.g. compared 
to the standard FRW cosmology with A as the additional phenomenological parameter. 

Perhaps the first to discuss light propagation in the clumpy universe was Zel'dovich 
in 1963 [1]. Before that, the idealized FRW geometry had been used as the basis in 
deriving theoretical predictions for the observable distance-redshift relations. Later in the 
go's, Zel'dovich's idea was elaborated at least by Bertotti [2], Gunn [3] and Kantowski [4]. 
However, it was the works by Dyer and Roeder in the 70's that made the idea more famous 
[5,6]; consequently, the distance-redshift relations allowing for the effects of clumping of 
matter on light are usually referred to as the Dyer-Roeder relations. 

It seems that the interest towards the effects of dumpiness on light propagation suffered 
a hiccup after 1976, when Weinberg argued against the use of the Dyer-Roeder version of 
the distance-redshift relations in place of the conventional FRW ones [67]. He proposed 
that the effect of light propagating in the empty intergalactic space would, on average, be 
canceled by an equal but opposite effect due to gravitational lensing caused by the matter 
clumps. However, Weinberg's argument has later been challenged [68]. In addition, we 
want to point out that the original Dyer-Roeder relations do not take two crucial physical 
properties of the universe into account: the inhomogeneities in the expansion rate and 
the growth of the inhomogeneities. Indeed, in Sect. ^ we demonstrate that the evolving 
inhomogeneous expansion can induce much stronger corrections to the observable distance- 
redshift relations than the mere dumpiness of matter in the original formulation of Dyer 
and Roeder - and hence unlikely to get counterbalanced by the gravitational lensing. 

2.5 The original Dyer-Roeder method 



Before presenting the generalized version of the Dyer-Roeder method in Sect. 3.2, we in 



troduce the original formulation here for reference and to justify the modifications. For 
more discussion on the method, see Refs. [5,6,69-72]. 

A good starting point for quantitative considerations of light propagation in the inho- 
mogeneous universe is the exact Sachs optical scalar equations [73]. They correspond to 
the scalar part of the Einstein tensor equation on null geodesies. The first one reads as 

f)f) 1 

^ + -e^ + a^ + R^uk^k'' = 0, (2.6) 

where 6 = V~^dV/dX is the expansion of the bundle of light rays taking up a volume 
V; A is the affine parameter along the null geodesies x^(A) that gives the energy of the 
photons as measured by an observer with four-velocity n^ via E{X) = —g^^u'^dx'^/dX; R^^ 
is the Ricci curvature tensor, whereas a^ represents optical shear and /c^ = dx'^/dX is the 
wave-vector of the light rays. The other Sachs equations are trivially satisfied when the 
optical shear is neglected as in the Dyer-Roeder approach and which is also an in-built 
assumption of the standard FRW description. 

To rewrite Eq. (2.6) in terms of observables, we use 6 = A~^dA/dX where A is the 



cross-sectional area of the light bundle and collect the derivatives into one term as 

to obtain 

^^^d^ + STTGT^.k'^k'' = , (2.8) 

where the Einstein equation has been employed in the form R^^, = SttGT^^ + g^^R/2 
with the facts that k^ is null, g^^pk^^k'^ = 0, and that the angular diameter distance is 
proportional to the square root of the cross sectional area of the light bundle, dA oc yA. 

By using ideal fiuid to describe the energy content of the universe, we have the energy 
tensor T^^ = {p+p)u^u^ +pg^^'^ where u^ = 5q in the comoving coordinates and p = Pm+Pr 
stands for the energy density of matter plus radiation and p = p^ is the radiation pressure. 
Hence, with the help of g^^ui^k" = -E{X), Eq. (|2^) yields 

^dAX) + 47TG{p + p)Ei{l + zfdAiX) = , (2.9) 

where we have used the definition of the redshift, E[X)/Eq = 1 + 2;, with Eq as the observed 
energy of the photons today. We still want the affine curve parameter A in terms of the 
non- affine but observable redshift z: 

where H{z) is the Hubble expansion along the path of the light and it has been assumed 
that the relation between the auxiliary scale factor and the redshift, a(0)/a(z) = 1 + 2;, 
is a good approximation and that the proper time t is independent of the spatial location 



(which is only possible if the vorticity of matter vanishes). Finally, by using Eq. ( 2.10 ) we 
can write Eq. ( |2.9| ) as 



d 
H{z)- 
oz 



{l + zfH{z)^^dA{z) 



47rG[p(z)+p(z)]dA(2)=0, 



(2.11) 



where the dependence on the energy of the photons Eq canceled out. 

By using the general relation between the angular and luminosity distances, 

dL{z) = {l + zfdA{z) , 



(2.12) 



which holds for geodesic light in any spacetime as proved by Etherington in 1933 [52, 74, 75], 
Eq. ( |2.11| ) can be transformed into the corresponding equation for the luminosity distance: 



{l + zfH{z) 



d_ 
dz 



[1 + zYH(z, 



d diiz) 



+ 4TTG[piz)+piz)]dLiz)=0 



(2.13) 



dz{l + zy_ 

The essential point Dyer and Roeder made was that, due to Eqs. (2.11) and ( 2.13| ) 
determining the optical properties of the universe, the matter density pm should be in- 
terpreted as the average density of the intergalactic regions where the observed light has 
traveled. To describe this, they introduced a constant smoothness parameter a, represent- 
ing the fraction of the matter density along our line of sight pm to the volume-averaged 
matter density pm of the whole universe in the FRW description. That is, pm = apm in 
Eqs. ( ^.11 ) and ( |2.13| ), with the value a = corresponding to a late universe where all 
the matter is concentrated into opaque clumps and the value a = 1 to the smooth early 
universe. Clearly, a can then be something between these two extremes. 

Eqs. ( p. 11 ) and ( p.l3|) are second order differential equations for dA{z) and diiz) with 



the initial conditions: dyi(O) = 0, ^^'(0) = Hq and ^1,(0) = 0, dL'{0) = H^ , where 
' = d/dz. Thus, we can easily write the solution to Eq. ( |2.13D as a power series in z: 



dT{z) 



H^ 



z + -z^ + 
4 



1 1 



a 



z^ + Oiz^) 



(2.14) 



where it has been assumed that the Hubble expansion has the flat FRW form with pres- 
sureless matter as the only source, H{z) = Hq{1 + z)^'"^, since the radiation is negligible 
at small redshifts where the solution ( |2.14| ) would only hold anyway. 

By comparing the Dyer-Roeder luminosity distance (|2.14| ) with its counterpart in the 
homogeneous ACDM model (^.5|), we see that there is no way a could mimic A, as it does 
not even appear in the all-important second-order term. In fact, the effect of clumping via 
a < 1 seems to effectively increase the deceleration of the expansion in concordance with 
the original results of Dyer and Roeder [5, 6]. 

However, on physical grounds, one could expect just the opposite result: due to the 
local gravitational attraction of matter, a shearless and irrotational region of space with 
low mass density should decelerate less than a denser one. This can also be seen from the 
Raychaudhuri equation, an exact scalar part of the Einstein equation [52]: 



1 d' 

yi/3 dT^ 



yi/3 



AttG 



-r 



^') 



(2.15) 
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where o"^ and uP' represent the shear and vorticity of matter respectively, r is the proper 
time of the matter particles and y is a local volume element, making V^''^ a generalized 
scale factor. As expected, in a region far enough from the high-density filaments so that 



UJ 



<C 'i'KGpm, the deceleration or the left hand side of Eq. (2.15) is the smaller the 



lower the local matter density pm- 

We argue that these physically reasonable conclusions do not show up in the solution 



(2.14) due to the incomplete description of the inhomogeneities in the original Dyer-Roeder 
method. Firstly, matter density pm is not the only quantity with inhomogeneous distribu- 
tion, but the local Hubble expansion along our line of sight H{z) can differ from the FRW 
value as well. In fact, the effect of the expansion rate on the distance-redshift relations 
is more dominant than the effect of the matter extent, because the expansion is mainly 
responsible for the redshift. This has been demonstrated within the LTB model in Ref. 
[56] and is manifest already in the distance-redshift relations of the FRW models: 



diiz) 



il + z) 



Ho' 



VW^^i 



--Sk 



V\l-^o\Ho 



dz 



(2.16) 



with only a relatively weak dependence on the overall density f^o via 



Sk{x) 



sinh(a;) 

X 

sin(x) 



if f7o < 1 
if f7o = 1 
if f7o > 1 



(2.17) 



whereas the expansion integrated over the past light cone, J^ H^^[z)dz, yields the main 
contribution to the relation ( 2.16| ). 

Altogether, we propose that in order to make the method physically more correct. 



the Hubble expansion H{z) in Eqs. ( 2.11 ) and ( 2.13 ) should be interpreted as the average 
expansion rate of the regions where the detectable light travels. For this, we introduce 
another parameter /?, the ratio of the expansion rate along our line of sight H{z) to the FRW 
value H{z), so that H{z) = j3H{z) in Eqs. ( ^.11 ) and ( 2.13 ). Another crucial ingredient 
missing from the original formulation is that structures evolve in the universe; that is, both 
of the parameters should depend on the redshift: a{z) and (3{z). The redshift-dependence 
of a was in fact already proposed by Linder in 1988 [69], but perhaps because the corrections 
in the original Dyer-Roeder method are so weak (see Eq. ( p. 14 )), the redshift dependence 
of a has won very little attention in the literature [76]. Instead, as we demonstrate in Sect. 
3.2, it is only after taking into account both the parameter (3 and its dependence on the 
redshift that the results start to go more hand in hand with intuition. 



3. Optically altered universe due to nonlinear structures 

In this section, we introduce the principal model of this work, where light propagation in 
the real clumpy universe is described phenomenologically by a generalized Dyer-Roeder 
method. Within this approach, the universe can still be homogeneous on large scales, 
but the nonlinear inhomogeneities on smaller scales render its optical properties different. 
As discussed in Sect. 0, the optical overall effect of the inhomogeneities can result from: 
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1) biased light propagation through voids that expand faster than the average, 2) global 
increase in the volume-averaged expansion due to spatial variations in the expansion rate, 
3) our possible special spatial location in the universe, or any combination of the three. 
There can also be effects that are not reducible to any of these; e.g. a clock-rate difference 
between underdense and overdense regions discussed in Refs. [40-45]. 

We proceed by summarizing the physical foundations of the model qualitatively in 



Sect. 3.1 and leave the discussion of its quantitative properties to Sects. 3.2 - 3.5. 



3.1 The physical foundations of the model 

The assumptions that form the physical foundations of the model, can be summarized in 
three main points: 

1. The older the universe, the clumpier it is and hence the emptier and larger the voids. 

2. The light from the well-observed objects travels mostly in the voids. 

3. The emptier and larger the voids, the faster they expand relative to the average. 

The first point is the standard notion of structure formation: initially overdense regions 
accrete more material via gravitational attraction and thus become denser and denser at 
the expense of the initially underdense regions, which thereby become emptier. Some 
observational evidence for the emptiness of voids has been found in galaxy surveys that 
show there to be typically < 1% galaxies near the center of the voids relative to the average 



galaxy density [7]. The second assumption was explained in Sect. 2^: the emptier a region 
the easier it is to see an object through it. The point three is an integrated effect of Eq. 
(2.15): denser regions have more mass and consequently stronger gravitational attraction 



to slow down the local expansion. 

The combination of these properties suggests a physical explanation for how the non- 
linear structure formation can be perceived as accelerated expansion without dark energy: 
as the detectable light traverses emptier and emptier regions, the expansion rate along the 
path of the photons increases relative to the average, causing accelerated expansion along 
our line of sight without local acceleration. More specifically, we present the following 
conjecture: due to the formation of voids and opaque structures, the average matter den- 
sity on the path of the light from the well-observed objects changes from Qm — 1 in the 
homogeneous early universe to Qm — in the clumpy late universe, so that the average 
expansion rate increases along our line of sight from EdS expansion Ht ~ 2/3 at high 
redshifts to free expansion Ht ~ 1 at low redshifts. 

Optically, the proposed phenomenon bears close resemblance to the local Hubble Bub- 



ble, discussed in Sect. |2^; see in particular Eqs. (|2.4D , (|2.5D and the paragraph thereafter. 
However, unlike in the models with a large void in the local matter distribution, it is our 
special location in time, not in space, that is here crucial for the effect, so the cosmological 
principle is not violated. On the contrary, since every observer in the late universe sees 
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the history of the structure formation, the phenomenon should look essentially the same 
no matter at which point the observer sits in the universe^. 

3.2 The generalized Dyer-Roeder method 



In Sect. 2.5, we reviewed a quantitative method to allow for the optical effects of the 
dumpiness of matter in the universe, following the reasoning of Dyer and Roeder [5, 6]. The 
basic idea of this method is to employ the average density along the paths from the objects 
to the observer in calculating the distance-redshift relations for those objects, instead of 
using an average density of the whole universe like in the FRW description. Albeit a 
promising starting point, we argued that the original method falls short of capturing two 
crucial physical phenomena: inhomogeneities in the expansion rate and the growth of the 
nonlinear structures. To take these properties into account, we present a generalized version 
of the original Dyer-Roeder method. 



Following the argumentation of Sect. 2.5, we take into account both the inhomo- 
geneities in the matter distribution and in the expansion rate via the redshift-dependent 
parameters a{z) and P{z). Continuing to neglect the vorticity of matter, we still have a 
coordinate system where the matter four- velocity is just u^ = 5q. Thus, we can simply 
substitute Pmiz) -^ a{z)pmiz) and H{z) -^ (3{z)H{z) in Eqs. ( |2.11| ) and ( p. 13 ) to obtain 



l5{z)H{z) 



d_ 

dz 



{l + zfl3{z)H{z)-^dA{z) 



{l+zfp{z)H{z) 



d_ 

dz 



dz 



{l + zfp{z)H{z) 



+ A-KG 



a{z)p„ 



4 

+ -^Pr 



z) 



dA{z)=f) (3.1) 



dz{l + zY 



+47rG 



a{z)pm{z) + -Pr{z) 



ddz) = , 

(3.2) 

where the symbols with an overbar represent the auxiliary FRW quantities and we assume 
that radiation does not cluster: Pr{z) = Pr{z). The initial conditions for the differential 
equations (|3]|) and (|]|) are: 4(0) = 0, 4'(0) = i?o~^ ^^^re x e {A,L]. 

To work out how the inhomogeneities can modify the luminosity distance in the gen- 



eralized Dyer-Roeder method, consider the solution to Eq. (3.2) as a power series in z: 



HodLiz) 



2/3(0) r + 



/3'^(0) P'{0) p"{0) a{0) 



3/32(0) 3/3(0) 6/3(0) 4/32(0) 



z-^+0(z^ 



(3.3) 



where radiation has again been neglected. Here the nonlinear structures, encapsulated in 
the functions a{z) and /3(z), have much more notable role than in the original Dyer-Roeder 
method in Eq. ( |2.14| ). In addition, as can be seen by comparing the solution ( p. 3D to the 
luminosity distance in the flat ACDM model ( |2.5D , the effect of the inhomogeneities is here 
importantly in the correct direction, mimicking A with /3'(0) < 0. 

We emphasize that the solutions dA{z) and diiz) to Eqs. (|3.lD and ( ^I^ ) do not repre- 
sent mathematically exact solutions of general relativity. Instead, their aim is to capture 
the most essential effects of the inhomogeneities on the observable properties of light that 



*This, at least, seems to distinguish the light propagation in freely expanding voids as the cause for the 
acceleration from Wiltshire's proposal whereby it is crucial that the observer sits in a static or gravitationally 
bound location, such as inside a galaxy, to perceive the inhomogeneity induced acceleration [40]. 
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the simplified FRW description fails to describe. The crudest approximations in Eqs. (| 



and (3.2) correspond to employing a different FRW result for each redshift. Even ap- 
proximations of this kind have proven to be accurate in calculating the distance-redshift 
relations [36]. Nevertheless, in a more thorough treatment, e.g. the utilized FRW relation 
a(0)/a(z) = 1 + z should be replaced with a corrected inhomogeneous result; see [25]. 

To calculate the observable distance-redshift relations dA{z) and ^^(z) from Eqs. ( pTl] ) 
and ( |3.2| ), we need the functions a{z) and P{z). An entirely phenomenological approach 
would be to determine them purely from the cosmological data analysis. In contrast, an 
ideal approach would be to derive a{z) and (3{z) from the full Einstein equations starting 
with some initial perturbations in the early universe. In this work, we employ an analytic 
approach which falls somewhere between these two extremes: based on theoretical and 
observational reasoning, we deduce a form for the functions a{z) and fi{z) containing a 
single phenomenological parameter that, given the approximations, would best correspond 



to the structure formation of the real universe; this is done in Sect. 3.3. In Sect. 3.4, we 
then fit the leftover parameter to various cosmological data. 

3.3 A phenomenological description of the nonUnear structure formation 

As boundary conditions given on a 3-dimensional hypersurface, the expansion rate H and 
the matter distribution pm are independent (apart from some constraints). In theory, one 
could thus give them almost any profile e.g. along our past light cone; that is, choose a{z) 
and f3{z) in Eqs. ( ^.1[) and (|3.2| ) arbitrarily. However, as suggested by the nearly uniform 
temperature of the CMB radiation, it is likely that the early universe has been very close 
to homogeneous. Thus, we only consider inhomogeneities that grow with time, which sets 
a constraint between the functions H{z) and Pm{z), or equivalently between a{z) and f3{z), 
reducing the number of independent functions to one. 

Some properties of a(z) and (3{z) can be deduced without resorting to a specific model. 
Firstly, according to Eq. ( |2.15| ), initial underdensities evolve into regions that eventually 
begin to expand faster than the global average. As the detectable light is assumed to 
propagate mostly in these regions, we have a < 1, hence implying that /3 > 1. The 
requirement of positive mass density sets a > 0, while Eq. ( 2.15| ) implies the upper bound 



P < 3/2. For an explicit constraint between a{z) and P{z), a specific model is required. 

The question for which we seek an answer is then: how much faster does a void with 
known matter density pm < Pm expand than the FRW value H^ i.e. given < a < 1, what 
is 1 < /3 < 3/2? We estimate this constraint by employing a different FRW result to each 
patch with different density and expansion rate. Given that we use the proper time of 
dust as our time coordinate, the constraint can be obtained by requiring a constant value 
for the age of the universe io between the different regions. This sets the matter density 
and expansion rate equal at t = for every patch, i.e. the condition that the universe was 
initially homogeneous. We do not, however, employ this method to model the dynamics of 
the universe as, e.g. due to vanishing vorticity and ignored mass flows, it cannot realistically 
describe the evolution of the structures at least when virialization becomes important. 

As a starting point, we take the functions a{z) and (3{z) to measure deviation from a 
flat FRW model with matter and radiation. Thus the quantities with an overbar in Eqs. 
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obey the relations 



47rGpm{z) = -Hl(lm{l + z) 



(3.4) 



H{z) = Ho{l + zf/^Vl + z^r , (3.5) 

where Clm = 87rGp,„(0)/(3^o) ^^'^ ^r = 87rGpr(0)/(3^Q) are the density parameters of 
the auxihary FRW model (with the constraint Ctm + ^r = 1) and Hq would be interpreted 
as the average global expansion rate in the FRW description. However, if the dynamical 
backreaction (|2.2| ) is significant, in the nonlinear late universe Pm{z) and H{z) are just 
auxiliary quantities with no obvious physical interpretation. 

To estimate the constraint between a{z) and (5{z)^ we first replace them by the auxiliary 
functions Hq{z) and (lm{z), defined through the equations: 



a{z) 



H^{z)nm{z) 
ff20 



(3.6) 



Piz) 



(3.7) 



With the definitions (|^ 

4:1TGpm{z) 



Ho{z) 1 + Z^rajz) + z{z + 2)nr 
Ho V {l + z){l + znr) 

), the line-of-sight quantities Pm{z) and H{z) take the form 



3 

4TTGa{z)pm{z) = -Hi{z)nn,{z){l + z] 



H{z) = f3{z)H{z) = Ho{z){l + z)Jl + zn^{z) + z{z + 2)i} 



{3.1 



(3.9) 



where, by construction, the deviation from homogeneity appears as the redshift dependence 
of Hq{z) and ^rn{z) which in the limit z ^ 1 reduce to the constants Hq and Vtm- Apart 
from this limit, there is no straightforward physical interpretation for Hq(z) and ilmiz)- 

The next step is to obtain the constraint between Hq{z) and Qm{z) by requiring that 
the age of the universe is the same between the different FRW patches. The relation 
between Hq{z) and ilmiz) then determines a constraint between a{z) and P{z) through 
Eqs. (|3.6D and ( ^.71 ). We are only applying the FRW results to voids, so 0^ < 1. In this 
case, the integration of the Friedman equation yields the age of the FRW dust universe: 



to 



\/l — ^m — J^marsinh 



l-Qr, 



(3.10) 



We use the redshift to label the different patches on the present-day spatial hypersurface, 
so Q,m and Hq in Eq. ( |3.10| ) become functions of z. The requirement for the age in Eq. 



( 3.10 ) to be equal between the different patches then yields Hq in terms of 0^ 



^o(^) 



1 — r2m(z) — r2m(z)arsinh 






(1 - f]™(z))3/2 






^^m\Z) 



(3.11) 
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where ig = constant is the age of the universe and in the last step we have employed an 
approximation which is accurate (|error| < 1.5%) in the interval < Q^niz) < 1- 

Using Eqs. ( |3.6D , ( |3.7D and ( |3.11| ), we can give 0.m{z) and Hq{z) in terms of a{z), 

/ I \ 2 

9 I . L 8 



n^{z) = - 1 - y 1 - -\lnraa{z) I (3.12) 



tn-^ 



H,{z) = ^ 1 + ^1 - -^(l„,a{z) , (3.13) 



so the line-of-sight quantities (p^) and (3.9) become 



inGp^iz) = ^t^2f^„a(z)(l + zf (3.14) 



^0^ 



H{z) = l3{z)H{z) = ^(1 + z) 1 + y 1 - - Vf^™«(^) 1 X 



\ 




1 + ^ l-\/l-^V^™«(^)l + z{z + 2){1 - n^) . (3.15) 



We impose two observationally motivated boundary values for the function a{z). 
Firstly, soon after decoupling the universe was still close to homogeneous, so the difference 
between the line-of-sight and globally averaged matter densities must have been negligible, 
corresponding to pm -^ Pm or a(llOO) ~ 1. At times earlier than z = 1100, the distinc- 
tion between the line-of-sight and volume averaged quantities does not make sense, since 
we cannot see beyond the opaque plasma wall formed by the last scattering surface. The 
second boundary value follows from the assumption that the regions where the detectable 
light from the well-observed objects travels today are nearly empty: Pm(0) <C Pm{0) or 
q;(0) ~ 0. Altogether, we then have a(0) ~ and a(llOO) ~ 1. 

With the endpoints q(0) ~ and q(1100) ~ 1 of the function a{z) fixed, we still need 
its behavior between these points. Once initiated, the structures can only grow in a universe 
with no dark energy, so a logical starting point is that a{z) must be a monotonically 
increasing function of redshift. We choose the exact behavior of the function based on 
physical understanding about the structure formation. 

A typical behavior of a small overdensity on a given scale L is that initially the excess 
decrease in its expansion rate is small, with linear proportionality to the initial condi- 
tions; but eventually, around a redshift zq{L), the gravity-driven cumulative increase in 
the clustering speed of the constituents collapses the system rapidly until it virializes due 
to pressure gradients and the conservation of angular momentum [77, 78] . 

A function that encapsulates these qualitative features is the exponential, a{z) = 
oq + ^^ Oje"^' ^*^^'' , summed over different scales. For simplicity, we consider only a single 
redshift zq when the voids start to dominate at the most significant scale. Taking into 
account the above deduced boundary values, we choose 

nmiz) = (l-e~'/'A^ , (3.16) 
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or 

a{z) = ^ (2 + e-"/^o) ^ (1 - e-"/^«) ^ . (3.17) 



Inserting Eq. (|3.17] ) in Eq. ( 3.13 ), we also have 



i7o(z)to^/i(^) = ^(2 + e-^/^«) . (3.18) 

In the end, we leave the transition redshift zq as a free parameter. A more sophisticated 
numerical estimation of the function a{z) will be performed in a future work, in which even 
{zi{Li)}i is a derived, not free, set of parameters. 

3.4 Comparing the model with cosmological observations 

In this section, we aim to find out how the model with a{z) given by Eq. ( 3.17| ) fits the 



main cosmological data sets: the cosmic microwave background anisotropy, the position of 
the baryon acoustic oscillation peak inferred from the galaxy distribution, the magnitude- 
redshift relations of type la supernovae, the local measurements of the Hubble constant 
and the Big Bang nucleosynthesis. For a thorough and open-minded review of the current 
observational status, see Ref. [79]. 

A good baseline for the studies is the model by Blanchard, Douspis, Rowan-Robinson 
and Sarkar [80,81], and especially its refined version by Hunt and Sarkar [82-84]. This is 
an FRW model with the following basic properties: no dark energy or A = 0, a mixture 
of cold and hot dark matter f^cDM = 0.8 & r^HDM = 0.1, flat spatial geometry /c = 0, 
a low Hubble constant h = 0.44 and a bump in the primordial perturbation spectrum 
at the scales k ~ 0.01. ..0.1 /iMpc~^, which could arise e.g. from phase transitions during 
inflation [85]; we shall call this the linear CHDM model. The model fits the CMB angular 
power spectra better than the standard ACDM model. It is also consistent with the 
nucleosynthesis constraints and with most of the features in the matter power spectrum of 
the galaxy distribution surveys. However, the model seems to fail in matching the observed 
magnitude-redshift relations of type la supernovae and the detected position of the baryon 
acoustic oscillation peak in the matter power spectrum and it is also in disagreement with 
the local measurements of the Hubble constant. For a more thorough description of the 
linear CHDM model, see [83]. 

A candidate for hot dark matter would be three species of massive neutrinos with 
Yli ''^i^i ~ 1-^ 6^' yielding Ohdm = ^u ^ 0.1. This is well below the present experimental 
upper bound of rriu = 2.3 eV [86], but is expected to be detectable in the forthcoming 
KATRIN /?-decay experiment [87]. Note that the sometimes quoted tighter upper bounds 
from cosmological observations depend heavily on the priors [88]. As neutrinos with mass 
higher than rriiy > 0.3 eV have been shown to cluster appreciably [89], their existence should 
not alter the assumption that the detectable light travels in almost empty space today. 

We consider a universe that is similar to the linear CHDM model at large redshifts, 
but instead of describing its optical properties with the standard FRW relations, we use 
the generalized Dyer-Roeder equations (|3.lD -( pl2| ) to determine the observable distance- 
redshift relations. Accordingly, the growth of small deviations from homogeneity is assumed 
to follow linear perturbation theory, so apart from the changes in the associated distance 
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measures, the perturbation spectra should be more or less unaltered. We call this the 
nonlinear CHDM model. We pose five observational tests for this model: 

1. To obtain the same value for the angular diameter distance to the last scattering 
surface as in the linear CHDM model. 

2. To match the position of the baryon oscillation peak in the matter power spectrum. 

3. To fit the observed magnitude-redshift relations of the type la supernovae. 

4. To explain the locally measured value of the Hubble constant. 

5. To remain within the limits on the baryon density set by the BB nucleosynthesis. 

By passing these tests successfully, the nonlinear CHDM model would explain the main 
cosmological data sets that support accelerated expansion. In Sects. 3.4.1| - pA6| , we discuss 



how the model survives these tests and the results are combined in Sect. 3.4.7 



3.4.1 The cosmic microwave background 

For the values D,b = 0.1, h = 0.44 and Oq = 1) the fitting formula of Hu and Sugiyama 
[90] yields Zdec = 1101.35 as the decoupling redshift, for which we use the approximation 
Zdec — 1100 throughout. 

The fit of the linear CHDM model to the WMAP 3-year data has been thoroughly 
performed by Hunt and Sarkar [83] so we only need to consider how the nonlinearities 
modify this. As the nonlinear structures appear only in late times, it is essentially only the 
comoving angular diameter distance dc{z) = (1 + z)dA{z) and the late integrated Sachs- 
Wolfe effects that can change. However, as the model already fits the CMB data without 
the nonlinear structures, it is desirable that adding the nonlinearities would not alter the 
distance to the last scattering surface, (ic( 1100). Physically, this is possible because the 
negative spatial curvature and the faster Hubble expansion in the voids have the opposite 
effects on the comoving angular diameter distance dc{z), so with a suitable value for the 
new parameter zq, their effects can cancel. 

By determining dc{z) numerically from Eq. (f3.1|), and dc{z) = dL^z) / [1 + z) from Eq. 



(2.16), we find that the value for zq that sets (ic(llOO) and dc{1100) identical, is zq = 0.347. 
Values in the range zq = 0.30. ..0.39 give less than 1% discrepancy between dc{lWO) and 
(ic(llOO). Although equal around z ~ 1100, the functions dc{z) and dc{z) are not equal at 
lower redshifts. 

Varying the Hubble constant does not change the limits on zq as Hq has essentially the 
same effect on both the "background" value (ic(llOO) and the value dc(llOO) which takes 
into account the nonlinearities. To conserve the successful angular power spectrum found 
in Ref. [83], we fix the FRW Hubble constant to the value ^o = 44 kms^^Mpc"^. 

Within the linear perturbation theory, a flat matter dominated FRW model does not 
have a late integrated Sachs- Wolfe (LISW) effect, because the perturbed gravitational 
potentials do not evolve in time, $ = 0. However, this does not necessarily hold in a 
universe with nonlinear structures if the CMB light has propagated mostly in voids. Indeed, 
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with $ potentially non-zero along our line of sight, one would expect a non-trivial LISW 
effect from nonlinear inhomogeneities also in a matter dominated universe with initially flat 
FRW geometry. On the other hand, as the LISW affects the angular power spectrum only 
at low multipoles /, its contribution to the goodness of the overall flt could be insignificant 
due to the large cosmic variance oc (2/ -|- 1)~^. Moreover, it might in fact be that the effects 
of the local structures dominate the angular power at the lowest multipoles [91, 92]. It is 
also possible that, even if the selection effects were significant in the supernova observations, 
they may be less important in the CMB observations. Altogether, along with other more 
involved issues, we postpone a quantitative discussion of the LISW to a future work. 

Apart from the potential modifications due to the LISW effect at small I, the angular 
power spectrum of the nonlinear CHDM model with zq = 0.347 is plotted in Fig. |l|. 
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Figure 1: Binned WMAP data points in the temperature power spectrum with theoretical 
predictions from the best fit CHDM and ACDM models; figure from Hunt and Sarkar [83]. 



3.4.2 Galaxy distribution 

As long as the nonlinear inhomogeneities only affect light, the evolution of the linear per- 
turbations remains the same as in the linear CHDM model. In the case of the CMB power 
spectrum this seems obvious, since the relevant evolution of the perturbations happened at 
z < 1100 when the nonlinearities were negligible. However, as the galaxy distribution also 
reflects the dynamics of the late universe, it is not clear that the nonlinear inhomogeneities 
would not play a dynamical role here. We nevertheless neglect the potential effects of the 
nonlinearities on the evolution of the perturbations, but remind that, in a more realistic 
treatment, these effects should be taken into account in the matter power spectrum. 
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As shown in Refs. [81,83], the hnear CHDM model fits the matter power spectrum 
apart from the position of the baryon acoustic oscillation peak. As in the CMB angular 
power spectrum, the main modification caused by the nonlinearities is the change in the 
distance-redshift relations. This opens up a possibility that the nonlinear CHDM model 
would fit the position of the baryon acoustic oscillation peak. To study this, we use the 
data from the Sloan Digital Sky Survey provided in Ref. [18]. 

The detected position of the baryon acoustic oscillation peak in the matter power 
spectrum constrains the distance to the redshift z = 0.35, corresponding to the acoustic 
oscillation scale today. More specifically, due to the different scaling of radial and angular 
distances, the exact distance measure constrained by the data reads as [18]: 

1/3 



Dy{z) = 



zdiiz) 



(3.19) 



H{z) 

As the actual observations are made in the redshift space and the distance-redshift rela- 
tion depends on the cosmological model, the "measured" values quoted in [18] are prior- 
dependent. We take this into account by using the more conservative 2a error limits. The 



value of the distance ( 3.19 ) to the redshift z = 0.35 deduced from the SDSS data is 

Dy(0.35) = 1370 ± 150 Mpc , (3.20) 

where the limits corresponding to 2a errors have been projected from Fig. 7 of Ref. [18] . 
With Hq = 66 kms~^Mpc~^ fixed, we find that the constraint set by Eq. ( p. 20 ) on the 



transition redshift parameter in Eq. ( |3.17| ) is zq > 0.25, whereas the best fit value from 



the CMB analysis of Sect. |X|, zq = 0.347, gives L>y(0.35) = 1480 Mpc. Fixing zq = 
0.347 requires Hq = 71 kms~^Mpc^^ or ^o = 47 kms~^Mpc~^ to give the best fit value 
1)^(0.35) = 1370 Mpc. 

In the case that the statistical properties of the primordial perturbations are the same 
at every spatial location in the universe, the baryon oscillation peak in the matter spectrum 
and the first peak of the CMB angular power spectrum should correspond to the same 
physical feature just seen at different spatial locations and at different epochs. The physical 
scale associated with these oscillations is set by the sound horizon at recombination. As the 
distance (|3.20 ) measures how large this length scale appears in the galaxy distribution, it is 



possible to combine it with the information of the associated scale of the peak in the CMB 

spectrum to make an additional observational test. The quantitative measure employed in 

Ref. [18] for this is defined as: 

_ IJy(0.35) 

^°-^^ = -djim ■ ^ ^ ^ 

It should be again noted that due to the prior-dependency of the distance-redshift 
relations, Eisenstein et. al. quote several "measured" values for -Ro.35- As dc(llOO) = 
13190 Mpc is the best fit value to the CMB data in the CHDM model, Eq. ( ^^ yields: 

i?o 35 = 0.104 ±0.011 , (3.22) 

where the limits should roughly correspond to 2a errors in the CHDM model. We find that 
the values of zq corresponding to these limits are: zq = 0.29... 0.83, with zq = 0.53 yielding 
the best fit value -Ro.35 = 0.104. As -Ro.35 is a ratio of distances, it is independent of Hq. 
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Albeit the nonlinear CHDM model with a{z) given by Eq. ( 3.17] ) is not as concordant 



with the SDSS data as with the other data sets, the value zq = 0.35 still falls within the 
experimental limits. Moreover, the results represent a huge improvement compared to the 
linear CHDM model, which yields the values Z)v(0.35) = 1760 Mpc and iio.35 = 0.134 that 
are both way off the rather conservative limits of Eqs. ( |3.2C| ) and ( |3.22| ). 

3.4.3 Type la supernovae 

We use the Riess et. al. gold sample of 182 type la supernovae [17] to study how the 
nonlinear CHDM model fits the SN observations. For this, we use the conventional x^ test: 

X' = ^( ^-^ ^^j , (3.23) 

n=l ^ " / 

where m{zn) = 5logiQ{dL{zn) /Mpc) + 25 is the theoretical prediction for the magnitude 
and an is the estimated error for the measured magnitude m (zn) of a source at the 
redshift Zn- 

Just like the fiat ACDM model with Qa and Hq as free parameters, the nonlinear 
CHDM model has two free parameters: the transition redshift zq and, say, the observable 
Hubble constant Hq. There is, however, a perfect degeneracy between Hq and the maximum 
intrinsic luminosity of the type la supernovae, Lja, so the uncertainty in the value of Lia 
refiects a similar uncertainty in the value of Hq constrained from this data. For this reason, 
the combination {Hq,Li^ is sometimes marginalized in the data analysis but we choose to 
instead use the fixed value for Lja employed by Riess et. al. in Ref. [17]. 

By letting both parameters vary, we find the best fit values Hq = 65.5 kms~^Mpc~^ 
and Zq = 0.39 that yield x^ = 161.1. A little inspection in the parameter space shows 
that there is degeneracy between zq and Hq, such that one can obtain good two-parametric 
fits (Ax^ < 5) roughly in the range: zq = 0.2. ..0.8. Keeping instead Hq fixed to the 
value 66 kms~^Mpc~^, we find that the best fit transition redshift is zq = 0.34, yielding 
X^ = 161.4. The increase Ax^ = 2.7 corresponds to 2a errors in a one-parametric fit 
and yields for the hmits: zq = 0.34±[]:J^. The absolute fit 161.4/182 = 0.89 is good and 
represents a tremendous improvement to the linear CHDM model where x^ = 1098 for 
h = 0.44 and even the best fit value h = 0.54 gives x^ = 283.4. 

For comparison, the best fit parameters of the ACDM model Oa = 0.67 and h = 0.63 
yield x^ = 158.8. However, for the concordance values Clx = 0.73 and h = 0.7 the fit is 
much worse: x^ = 305.2, but this can be fixed by changing the assumed value of Lia or 
equivalently changing h. Indeed, we find that the value h = 0.642 gives the best fit for 
the concordance Cl\ = 0.73 ACDM model with x^ = 162.5. Overall, the conclusion from 
comparing the results between the nonlinear CHDM model and the ACDM model is that 
the goodness of their fits to the Riess et. al. data are essentially equal. In Fig. |2|, we plot the 
residual magnitude-redshift relations for three different models: the linear CHDM model 
with Ji = 0.54; the ACDM model with ^a = 0.67 and h = 0.63; and the nonhnear CHDM 
model with zq = 0.34 and h = 0.66. 
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Figure 2: 16 binned data points from the Riess et. al. gold sample in the residual magnitude- 
redshift diagram with theoretical predictions from three different models: the linear CHDM model 
with h = 0.54 (dashed black curve), the ACDM model with CIa = 0.67 and h = 0.63 (dashed dotted 
blue curve), and the nonlinear CHDM model with zo = 0.34 and h — 0.66 (solid red curve). We 
have subtracted the h = 0.66 freely expanding Milne universe relation from the curves to produce 
the residual diagram, so only the h = 0.66 nonlinear CHDM model intersects the origin. 



3.4.4 The Hubble Key Project 



The best fit values obtained from the different data sets of Sects. 3.4.1 , 3.4.2 and 3.4. 3| for 
the observable Hubble constant are all consistent with the low redshift measurements of 
the expansion rate that should specifically probe the observable Hubble constant; the value 
Hq = 66 kms~^Mpc~^ sits firmly within the limits of the final results from the Hubble Key 
Project: i^o = 72 ± 8 kms^^Mpc^^ [93] or Hq = 62.3 ± 5.2 kms^^Mpc^^ [94]. In addition, 
there are persistent indications that one gets systematically lower values for the Hubble 
constant by observing objects at higher redshifts [79]. For the standard FRW cosmology 
this feature is an anomaly while in the nonlinear CHDM model it is a prediction: towards 
higher redshifts, the optically perceived expansion should be closer to the (lower) FRW 
expansion rate due to the fact that the earlier universe was more homogeneous. 



3.4.5 Age of the high redshift quasar APM 08279 + 5255 

The particularly well-observed high redshift quasar APM 08279 + 5255 [95, 96] at z = 3.91 
with estimated age 2... 3 Gyr has been proposed as a litmus test for every cosmological 
model [97]. For several cosmological models, most notably the homogeneous standard 
ACDM model where tage(3.91) = 1.7 Gyr and almost all the dark energy models, the 
observed quasar implies an age problem: the universe appears to be younger than its 
contents. 
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The considered nonlinear model offers a resolution to the age problem: if the observed 
light from the quasar propagates in voids that expand faster than the average, the light 
redshifts more rapidly than expected from the stretching of the global FRW scale factor, 
implying that the quasar is seen in an older universe than predicted by the age-redshift 
relation of the homogeneous FRW model. An order of magnitude estimate can be obtained 
by approximating that the space along our line of sight to the quasar has expanded freely. 
This increases the age at a given redshift z relative to the FRW value by a factor \/l + z, 
yielding tagc = 3 Gyr for the age of the universe at z = 3.91, consistent with the estimated 
age 2.. .3 Gyr of the quasar APM 08279 + 5255. 

3.4.6 The Big Bang Nucleosynthesis 

Finally, we comment on the nucleosynthesis. The overall result from deducing the abun- 
dances of the light elements from the spectra of stars and galaxies is that the baryon 
density falls within the range 0.017 < VIb^ < 0.024 at 95% confidence level [98]. As long 
as the early universe was nearly homogeneous and the constants of nature do not vary, 
these limits are model-independent. Therefore, the observed value should not be affected 
by the nonlinear inhomogeneities in the late-time universe so that the value for the baryon 
density, Clsh? = 0.1-0.44^ = 0.0194, falls within the above quoted observational limits also 
in the nonlinear CHDM model. 

3.4.7 A concordant nonlinear CHDM model with no dark energy 

By combining the results from the analysis of the different data sets in Sects, p. 4. 1 -3. 4.6, 



we obtain a concordant cosmological model with the following properties: the observable 
Hubble constant h = 0.66, the transition redshift to the era when nonlinear structures 
become dominant zq = 0.35, baryon proportion CIb = 0.1, cold dark matter proportion 
^CDM = 0.8, hot dark matter proportion ^hdm = 0.1 and no dark energy A = 0. All the 
data sets are consistent with these values and, apart from the SDSS data, the concordance 
value for zq is actually very close to the best fit value of each separate data set; see Table 
1. Moreover, the concordance value for Hq yields to = Hq = 14.8 Gyr as the age of the 
universe, consistent with the lower bound determined from the astronomical estimates for 
the age of the oldest globular clusters: 11.2 Gyr (at 95% confidence level) [99]. The main 
properties of the best fit nonlinear CHDM model are summarized in Table 2. 

Description CMB BAO peak Type la SNe Concordance value 

Transition redshift zp 0.35ifci]:[i| 0.53ifc[]:i^ O-^^^om 0-35 

Table 1: The constraints on the transition redshift from each data set and the concordance value. 

3.5 Discussion of the fits 

For consistency, we have tested various forms for the structure formation function a{z), in 
addition to the form ( ^.17 ) employed in the data analysis in Sect. 3^. The generic outcome 



seems to be that the exact form of the function a{z) is not crucial to fit the data as long 
as it changes roughly from q(1100) ~ 1 to a(0) ~ 0. Neither does the fit require that the 
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Description 


Parameter 


Concordance value 


Transition redshift 


zq 


0.35 


Observable Hubble constant 


h 


0.66 


FRW Hubble constant 


h 


0.44 


Age of the universe 


to 


14.8 Gyr 


Baryon proportion 


^B 


10% 


Cold dark matter proportion 


^CDM 


80% 


Hot dark matter proportion 


S^HDM 


10% 



Table 2: Best fit parameters for the nonhnear CHDM model with no dark energy. Only one of the 
three parameters h, h and io is independent. It is noteworthy that a phenomenologically similar 
model was found by Wiltshire et. al. [42], though with different interpretation for how the voids 
give rise to the acceleration [40-45] . 



observed light should travel in totally empty space today: it is possible to have nonzero 
present-day average matter density on our line of sight a(0) ~ O(O.l) and still fit the data 
sets, although with slightly different best fit values for zq, to and the new parameter a(0). 
An important implication is that the seemingly worrisome feature of the matter density 
along our line of sight becoming exactly zero today {<=^ a(0) = 0) is not necessary, but 
we have chosen the one-parametric form ( 3. IT] ) for simplicity. 
As an example, the more sophisticated two-parametric form 



a{z) 



1 -|- tanh(2;/zo 



5 — tanh(z/2;o — (p) 



(3.24) 



with the values zq = 0.35 and (j) = 1.46 yields an even more concordant fit to the combined 
data sets; it gives Dy(0.35) = 1370 Mpc for (l3!2C| ), i?o.35 = 0.104 for (|3^ , x^ = 163 for 
the goodness of the supernova fit and dc (1 100) /dc (1100) = 1.000 or no deviation from the 
distance to the LSS in the best fit linear CHDM model. 



Due to the approximations employed in Sects. 3.2 and 3.S to enable an analytic treat- 
ment, we expect that the numerical values of the best fit parameters can represent the effect 
of the inhomogeneities only at a relatively crude level. In a more realistic treatment, the 
analytic approximations should be replaced with more thorough numerical procedures and 
perform the fits to all the data sets from scratch - e.g. by making the necessary data cali- 
brations using a nonlinear model as the baseline instead of the linear ACDM model. This 
concerns especially the matter power spectrum, for which the full analysis requires consid- 
erably more effort and is therefore postponed to a future work. Instead, due to the mutual 
degeneracy between the nonlinear structure formation and the cosmological constant JIai 

the fit does not gain a noteworthy improvement by 



evident e.g. in Eqs. (|2.5|) and (|3^ 
keeping Q.\ in the analysis. 

An inevitable consequence of considering observables that represent averages over dif- 
ferent directions on the sky is the mathematical resemblance between the model of Sect. 
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m and a geocentric LTB bubble. Naturally, with photons from different directions prop- 
agating through different void profiles, the isotropic Dyer-Roeder distance-redshift rela- 
tions can only represent idealized ensemble averages. The actual observations should thus 
contain angular variations that are relatively smaller when making observations over dis- 
tances much greater than a homogeneity scale and hence the largest at low redshifts. In- 
deed, considerable amount of scatter with unidentified cause has been observed e.g. in the 
magnitude-redshift relations of the type la supernovae and importantly, especially at low 
redshifts [100-102]; for detailed analysis of scatter in HKP data, see [103]. The observed 
anomalies at the low multipoles of the angular power spectrum represent corresponding 
features in the CMB data [91,92]. Moreover, Gurzadyan et. al. have identified features in 
the randomness of the raw CMB data that could originate from voids [104]. Altogether, 
it appears a promising option that these apparent peculiarities would in fact be caused, 
at least partially, by the nonlinear inhomogeneities and thus provide evidence for their 
cosmological role addressed in this work. However, due to our lack of knowledge from the 
detailed (dark) matter distribution and expansion rate, it seems difficult to make precise 
quantitative predictions for the angular deviations. An important task is to search for 
correlations between the apparent anomalies of the different data sets. 

4. Conclusions 

We have studied the effect of structure formation on the observable distance-redshift rela- 
tions by generalizing the Dyer-Roeder method to allow for two important physical prop- 
erties of the universe: inhomogeneities in the expansion rate and the growth of nonlinear 
structures. We noted that light propagation through voids can cause the expansion along 
our line of sight to increase from flat dust FRW-expansion at high redshifts to free expansion 
at low redshifts and thereby have a similar effect on the observed distance-redshift relations 
as dark energy. Conjecturing such a transition, we showed that the inhomogeneities offer 
a natural physical explanation for the cosmological observations without dark energy even 
with a single physically well-grounded phenomenological parameter, the transition redshift 
zq to the era when the matter density becomes negligible along our line of sight or when 
the detectable light starts to propagate in freely expanding space. 



Indeed, applying the generalized Dyer-Roeder method in Sect. |3.4| , we found a simple 
inhomogeneous model that fits the observations from the CMB anisotropy, the position of 
the baryon acoustic oscillation peak inferred from the SDSS data, the magnitude-redshift 
relations of type la supernovae, the local Hubble flow and the Big Bang nucleosynthesis, 
plus is potentially consistent with the estimated age of the quasar APM 08279 + 5255 at 
z = 3.91. The main features of the resulting concordant model are: 90% dark matter, 10% 
baryons, zq = 0.35 as the transition redshift to the void-dominated era, the values h = 0.66 
and h = 0.44 for the observable and FRW Hubble constants that yield to = 14.8 Gyr as the 
age of the universe, consistent with astronomical estimates, and no dark energy of any kind 
or A = 0. At large redshifts, this nonlinear CHDM model is similar to the linear CHDM 
model by Hunt and Sarkar [82,83]. A non-trivial point in the best fit model is that the 
distance-redshift relation is not a reproduction of the corresponding relation in the ACDM 
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model. The model is based on inhomogeneities of the observed kind so, unlike a large void in 
the local matter density, there is no violation of the cosmological principle: it is our special 
location in time, not in space, that is required for the acceleration. Therefore, a prediction is 
that every observer in the present-day universe should perceive the inhomogeneity induced 
acceleration. 

Although homogeneous models can also provide a phenomenological fit to the cosmo- 
logical data by assuming dark energy, they entail various issues that can be formulated 
into five questions: 

1. What justifies the use of the FRW metric to model light propagation in the real 
universe? 

2. Why would dark energy appear at a very low temperature T ~ 4 K if it has had the 
huge temperature range from Tp ~ 10^^ K to To ~ 2.7 K available? An expectation 
for the probability of a quantum field to arise at a temperature T is P ~ e""^^/^. 

3. Why would the effects of dark energy appear when large voids start to form? 

4. Why would A be fine-tuned such that the universe undergoes nearly free expansion 
today, H()t() = J^ arsinhw -^_^ ~ 1, whereas a slightly bigger A would yield 
HqIq > 1? This issue was pointed out by Kutschera and Dyrda in Ref. [105]. 

5. If nonzero, why would the cosmological constant be so tiny, A ~ 10^"*^^^ G~^ ■, whereas 
naive dimensional analysis suggests A ~ G~^? 

The structure formation offers a smart solution to the coincidence problem: if induced by 
the voids, the onset of the perceived acceleration naturally coincides with the formation of 
the voids. Furthermore, the observed -ffo^o ~ 1 corresponds to the characteristic expansion 
rate of freely expanding voids, t ~ to/2 is a natural time scale for gravity to form structures 
on the relevant length scales and A = is compatible with the data if void formation is 
responsible for the acceleration. In addition, there are various observed features in the 
universe that are anomalies for the standard FRW cosmology, such as the cold spot in the 
CMB, the low multipoles of the CMB, the CMB ellipticity and the scatter in the supernova 
data, but could be explained naturally by the anisotropy of the cosmic void distribution 
in the sky; a more exhaustive list of such features can be found in Ref. [44]. As we have 
only discussed how to construct models of the observed inhomogeneous universe within 
relativistic classical gravity, the answer to the long-standing debate of why would A = is 
naturally beyond the scope here. 

The transition to the nonlinear era, roughly given by zq = 0.35, might seem to happen 
rather late, only At = 5.4 Gyr back in time, whereas the oldest stars have formed already 
At > 12 Gyr ago. However, there are in fact several natural reasons to expect such a seem- 
ingly late transition. Firstly, while traveling through structures that are still collapsing 
in earlier times, light is blueshifted which can partially counterbalance the excess redshift 
caused by the light propagation in voids. Later, however, when the structures have viri- 
alized also at larger scales, there are more and more structures which neither collapse nor 
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expand. Those photons that are able to penetrate these viriaUzed overdense regions, then 
do not blueshift anymore, diminishing this partial counterbalance at late times. Secondly, 
the large amount of weakly interacting dark matter in the universe prefers to form large 
halos which cluster into smaller clumps later than the luminous matter which nonetheless 
indicates the clumping to us. It also appears likely that the voids can grow large enough 
only in the late universe. Finally, when contrasting the acceleration from structure for- 
mation with dark energy driven acceleration, it is essential to note that the rate at which 
gravitational phenomena happen is largely independent of the temperature, whereas non- 
gravitational interactions happen faster at higher temperatures. Therefore, while it seems 
extremely difficult to explain why a quantum field would arise at T ~ 4 K especially if 
it has had the huge temperature range from Tp ~ 10^^ K to Tq ~ 2.7 K available, it is 
natural that the inhomogeneity induced acceleration would start as "late" as At ~ 7 Gyr 
ago (corresponding to redshift z ~ 0.52 that can be read off e.g. from Fig. |2|). 

The analytic approach employed in this work necessitates approximations, so the pre- 
cision of the numerical values for the best-fit parameters presented here is not expected to 
be as high as an ideal analysis could yield. Thus, it is difficult to say anything conclusive 
about the subtler features of the concordant nonlinear CHDM model, such as the bump 
in the primordial spectrum or the small hot dark matter component. Instead, given the 
physically well-grounded foundations of the proposed conjecture about light propagating 
mostly in voids, the observational evidence for the nonlinear structures as well as the large 
volume taken up by voids, the concordant fit to the cosmological data and the offered smart 
solutions to some of the perhaps most puzzling theoretical problems in the modern cos- 
mology, it seems perfectly conceivable that the main physical implications - most notably 
that dark energy may not exist - could be correct. 

An important future test is to study quantitatively the observed angular deviations in 
various data that are an anomaly for the homogeneous standard ACDM model, but a pre- 
dicted feature in the inhomogeneous models arising from the anisotropic void distribution. 
Another crucial test is to demonstrate the proposed conjecture within less approximative 
solutions of general relativity. Ideally, a thorough derivation of the conjecture from theory 
would involve a numerical computation employing the full Einstein equations starting with 
some initial perturbations in the early universe. More realistically, the feasible next step is 
to replace some of the analytic approximations with more careful numerical procedures; we 
leave this analysis to a future work, where e.g. the transition redshift zq will be a derived, 
not phenomenological, parameter. Although potentially laborious, these tests should be 
realizable with present-day knowledge, so the proposed conjecture is certainly falsifiable. 
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